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Abstract 

We analyse W boson production at hadron colliders in association with one or 
two jets, both with the exact kinematics and in the high-energy limit. We argue 
that the configurations that are kinematically favoured tend to have the W boson 
forward in rapidity. Thus W boson production in association with jets lends itself 
naturally to extensions to the high-energy limit, which we examine both at leading 
order and by resumming higher-order corrections through the BFKL theory. 



1 Introduction 



In recent years strong-interaction processes characterised by two large and disparate en- 
ergy scales, which are typically the squared parton center-of-mass energy s and the squared 
momentum transfer t, with s ^> i, have been analysed. These processes can be divided 
into two categories: a) inclusive processes, such as dijet production in hadron collisions at 
large rapidity intervals [l], forward jet production in DIS ||, |], ||, and 7*7* collisions in 
double-tag events, e + e~ — > e + e~+ hadrons ||; b) diffractive processes, such as dijet pro- 
duction with a rapidity gap between the tagging jets, either in hadron collisions [[/], H, || 
or in photoproduction |L(| . 



The interest in these processes stems from the possibility that their description in terms 
of perturbative-QCD calculations at a fixed order in the coupling constant a s might not 
be adequate, and that a resummation to all orders of a s of large contributions of the 
type of ln(s/|t|), performed through the BFKL equation [IT, I2~|, might be needed. An 
additional motivation for the analysis of processes in the limit s>t, and in particular for 
dijet production in hadron collisions at large rapidity intervals, inclusive or with a rapidity 
gap, is to use it as a test ground for the production of a Higgs boson in association with 
jets at the LHC. A Higgs boson is mainly produced via gluon fusion, g g — > H, mediated 
by a top-quark loop. If the Higgs-boson mass is above the threshold for vector-boson 
production, the Higgs boson decays mostly into a pair of W or Z bosons. The signal, 
though, is likely to be swamped by the WW, QCD and tt backgrounds. A Higgs boson 
of such a mass is also produced in q q —>■ qq H via electroweak-boson fusion, W W and 
Z Z — > H, though at a smaller rate [|13| . However, this would have a distinctive radiation 
pattern with a large gap in parton production in the central rapidity region, because the 
outgoing quarks give rise to forward jets in opposite directions [jy], [151 , with no colour 



exchanged between the parent quarks that emit the weak bosons |16| 17] . Accordingly, 



the topology of the final state has been used to reduce the overwhelming WW + 2-jet 
background fL8| . In fact, requiring in W W + 2-jet production two forward jets in opposite 
directions, which entails a large dijet invariant mass, makes the parton sub-processes to 
be dominated by gluon exchange in the crossed channel, with the W's produced forward 
in rapidity. 

In this paper, we analyse forward W production in association with jets as a natural 
extension of dijet production at hadron colliders and forward-jet production in DIS, and 
as a process that for large dijet invariant masses shares the same dynamical features 
(i.e., gluon exchange in the crossed channel) as WW + 2-jet production with forward 
jets, but is considerably simpler to analyse. There are additional reasons to consider 
this process: firstly, it could be experimentally easier to pick up forward W bosons that 
decay leptonically than forward jets; once a forward lepton has triggered the event, one 
observes the jets that are associated to it, with no limitations on their transverse energy. 
Conversely, in a pure jet sample one usually triggers the event on a jet of relatively high 
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transverse energy, thus the triggering jet cannot be too forward. Secondly, W production 
in association with jets lends itself naturally to extensions to the high-energy limit, since 
it favours configurations with a forward W boson, as we shall see in Section [2.2| . We 
limit our analysis to H^-boson production, however we expect the same kinematical and 
dynamical considerations we make in this work to apply to Z-boson production as well. 

In Section Q we examine the exact leading-order inclusive rapidity distributions for 
W + 1-jet and W + 2-jet production, for each parton subprocess, as well as the rapidity 
distribution of the W boson when the rapidity interval between the two jets is large. 
In Section ^ we review the high-energy factorisation and the derivation of the impact 
factor for jet production, and we calculate the impact factor for W + 1-jet production. 
In Section |] we consider the rate for W + 2-jet production in several high-energy ap- 
proximations, and in Section [5] we review the BFKL Monte Carlo event generator. In 
Section |] we analyse several distributions and candidate BFKL observables, such as the 
production rate as a function of the rapidity interval between the two jets, the azimuthal 
angle decorrelation and the mean number of jets. Finally, in the last section we draw our 
conclusions. 



2 Kinematics of W + 1-jet and W + 2-jet Production 

In this section we analyse in detail the kinematics of W production in association with one 
or two jets, and we show that in pp colliders asymmetric configurations with a forward 
W boson are naturally favoured.^ The results presented here have been obtained using 
tree-level matrix elements generated by MADGRAPH ]19[ . 



2.1 W + 1-jet production 

We consider the hadroproduction of a W boson with an associated jet. At leading order 
(LO), the parton subprocesses are qq — > W g and qg — ► Wq. The momentum fractions 
of the incoming partons are given through energy-momentum conservation by 

_ \Pj± I Vi , m ^ r y w 
V s V s 

Xb = IPiAe-Vi + ^e-vw^ ( 2.i) 



with pj x the jet (and the W) transverse momentum and m± = ym^ + |p J± | 2 the W 
transverse mass. 

What are the typical distributions in yj and in y w l At proton-antiproton colliders, 
the subprocess q q — > W g is leading; since the incoming quark and antiquark are valence 

* Unless stated otherwise, we always understand W to include both W + and W~ production. 
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quarks and the up and the down quark distribution functions have different shapes, this 
entails an asymmetry in the rapidity distribution of W + versus W~ bosons, both in 
fully inclusive (Drell-Yan) W boson production [^U] and in W + 1-jet production, and 
accordingly a large plateau for the rapidity distribution of the W boson as a whole. 

Also at proton-proton colliders, the W boson may be produced abundantly in the 
forward rapidity region. As in the rapidity asymmetry, the physical mechanism is 
the difference in the shape of the p.d.f.'s of the incoming partons. In fact, to be definite 
let us consider the subprocess q g — > W q, which at proton-proton colliders is dominant, 
and suppose that the incoming gluon enters from the negative-rapidity direction while the 
quark enters from the positive-rapidity direction, so we can identify x a as the gluon and 
Xb as the quark momentum fractions. The gluon distribution function is very steep, so 
it pays off to have x a as small as possible, which can be achieved by taking y w negative. 
That increases considerably the value of x&, however, because of the shape of the valence- 
quark distribution function, it can be achieved without paying a high price. In Fig. [I], 
we consider the rapidity distribution of the W boson in W + 1-jet production, broken 
up into its parton components. The renormalisation and factorisation scales, fi R and fj, F , 
are taken to be equal to (\pj ± \ + m±)/2. In Figs. we have taken the W mass to be 
m w = 80.44 GeV, we have used the p.d.f.'s of the package MRST99cg and evolved a s 
accordingly | 2l| . Applying the argument above to both the gluon incoming directions for 
qg — > W q, yields a broad rapidity distribution of the W boson, Fig. |I]a. The picture 
above applies to the subprocess qq^W g too, Fig. [TJt> , since the antiquark is in this case 
a sea quark. 

In case the W boson is produced forward in rapidity, with which rapidity is the jet 
typically produced ? If the jet is produced in the opposite hemisphere with respect to the 
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Figure 2: Rapidity distributions of the jet for the subprocesses (a) qg — > W q and (6) 
qq W g at the LHC centre-of-mass energy = 14 TeV and with Pj ±m i n = 30 GeV. 



boson, the rapidity interval \y w — yj\ is large, however we know that in this case W + l- 
jet production is strongly suppressed (at LO), since its parton subprocesses can only have 
quark exchange in the crossed channel, and thus the related production rate falls off with 
the parton centre-of-mass energy s. Thus this configuration is dynamically disfavoured. 
On the other hand, jet production in the same hemisphere as the W boson or centrally 
in rapidity keeps x a small without substantially increasing x&. However, whether the jet 
is produced in the central region or in the same hemisphere as the W boson depends on 
the detailed shape of the p.d.f.'s, namely on how large we can afford to make Xf, while 
keeping x a small. In Fig. ^| we plot the rapidity distributions of the jet at LHC energies. 



2.2 W + 2-jet production 

Let us consider the hadropro duct ion of a W boson with two associated jets. At LO the 
parton subprocesses are 

(a) gg^Wqq, 

(b) qq-^Wgg + Wqq, 

(c) qq — > W qq , 

(d) qg^Wqg. (2.2) 

The momentum fractions of the incoming partons are given through energy-momentum 
conservation by 



x a = fa^ 1 + %^ + ^e 
V s V s V s 
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subprocesses 






gg -> Wqq 


170 


170 


qq^Wgg + Wqq 


580 


400 


qq — > Wgg 


400 


300 


g# -> Wqg 


3300 


2200 



Table 1: Total cross sections (pb) for the production of W ± boson in association with 
two jets with transverse momentum Pj 12± > 30 GeV and interjet distance R(ji,j2) = 
V (Vji ~~ l/j'2) 2 + {4>h — 0j 2 ) 2 — 0-4 on the rapidity- azimuthal angle plane. 



Xb = \EhA e -Vn + \PjjA e -yj 2 + O^e-vw , (2.3) 
\/s V s V s 

with pj 1 2± the jet transverse momenta and m± = \Jm\, + \pj 1± + Pj 2± \ 2 the W transverse 
mass. For the four subprocesses of Eq. ( f2.2| ), the total cross section for the production of 
a W boson in association with two jets is given in Table |l|. 

What are the typical rapidity distributions of the W boson and of the two jets? In 
Figs. ^|-|] we plot the rapidity distributions of the W and of the two jets in W + 2-jet 
production. The renormalisation and factorisation scales, fi R and /i F , are taken to be 
equal to (|pj 1± | + \Pj 2± \ + m±)/2. The subprocess gg — > Wqq is perfectly symmetric, 
thus the W boson and the two jets are produced mostly in the central rapidity region. 
However, in the other subprocesses that is not the case: looking at the distributions in 
y w (Fig. ^) we see that as we move from (a) to (d) the W boson tends to be produced 
more and more forward in rapidity. Examining the distributions in t/j 2 (Fig. [|), where j% 
is the jet that is closest to the W, we see that this jet tends to follow the W in rapidity. 
From the distributions in — yj 2 (Fig. ||), we see that in (a) and (6) jet 1 tends to be 
produced more centrally; in (d) it follows the W boson and jet 2, thus emphasizing the 
kinematical features already noted in W + 1-jet production (the twin peaks observed in 
Fig. I in (a), (b) and (d) are due to requiring two jets with interjet distance R(j\, J2) > 
0.4); finally in (c) it tends to be produced far in rapidity from the W boson and jet 2. 

To understand how these configurations come about, we consider qg Wqg and 
follow the analysis of Section [2"Tl"| , i.e. we identify x a as the gluon and Xb as the quark 
momentum fractions. To make x a as small as possible at the price of increasing Xb, the 



W boson is produced forward (Fig. |3|). Note that with respect to Section |2.1| this is made 
easier by the presence of two jets, which let the W boson have a transverse momentum 
as small as kinematically possible: ultimately, when the jets are balanced in transverse 
momentum, the W transverse mass reduces to the mass, m± — > m w . In addition, one 
jet, say j 2 , is always linked to the W boson via a quark propagator as in W + 1-jet 
production, so it tends to follow the W in rapidity, as in Fig. f|, however the position 
of the other jet is a dynamical feature peculiar of W + 2-jet production: thanks to the 
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gluon exchanged in the crossed channel, that jet can be easily separated in rapidity from 
the W boson. In qq — > W g g + W qq, the kinematical mechanism is the same as in 
qg — > W q g since the antiquark has a sea quark p.d.f., however only qq — > W qq can have 
a gluon exchanged in the crossed channel. For g g — > W qq, which has equal p.d.f. 's for 
the incoming particles and no gluon exchanged in the crossed channel, we obtain a central 
distribution, as expected. Note, however, that in Fig. | and following the contribution of 
g g ^ W q q to W + 2-jet production is quite small. The qq — > W qq channel is peculiar, 
since the largest contribution comes from valence-quark distributions, which tend to have 
rather large x's. In addition, at the dynamical level it features only diagrams with gluon 
exchange in the crossed channel. Thus to make one x large, it tends to have the W boson 
and a jet slightly forward in rapidity, while to make the other x large, it has the second 
jet well forward (and opposite) in rapidity. 

Next, we require that the two jets are produced with a sizeable rapidity interval, 
\Vh ~ Vh\ — 2, and look at the rapidity distribution of the W boson with respect to 
the jet average, y w — (y^ + yj 2 )/2 (Fig. |6p. Now the requirement that the rapidity 
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Figure 4: Distributions in yj 2 , where j 2 is the jet that is closest to the W, for the subpro- 
cesses of Eq. ( J2.2|) at the LHC centre-of-mass energy = 14 TeV and with Pj ±m m = 30 
GeV. 



interval between the jets is large makes the subprocesses with gluon exchange in the 
crossed channel stand out even more, and the W boson, which is linked to one of the 
jets by quark exchange in the crossed channel, to follow that jet in rapidity. This is 
stressed by the double peaks in (b), (c) and (d). Note that the dip between the peaks is 
maximal for qq — > W qq, which features only diagrams with gluon exchange in the crossed 
channel. Conversely g g — > W qq yields the W boson in the central rapidity region and 
approximately equidistant from the two jets, however it is strongly suppressed since it 
can only have quark exchange in the crossed channel. 

The plots of Fig. |6] are characterised by the dominance of the subprocesses featuring 
gluon exchange in the crossed channel. The same feature is exhibited by events where 
we select the jet, say ji, that in rapidity is furthest away from the W, require that 
\yw ~ Vji\ ^ 2, and examine the distribution in y w — yj 2 . Since in this case j2 is always 
linked to the W boson by quark exchange, the distributions are all centered about zero. 
The leading subprocesses factorise then naturally into two scattering centres: an impact 
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factor for W + 1-jet production, and an impact factor for jet production. The two impact 
factors are connected by the gluon exchanged in the crossed channel. Accordingly, the 
dashed lines of Fig. || have been obtained by taking the high-energy limit of the amplitudes 
featuring crossed-channel gluon exchange (see Eq. (|4.3|) ). On these amplitudes we can then 
insert the universal leading-logarithmic corrections of 0{a n a ln n (s/|t|)), and resum them 
through the BFKL equation. The impact factor for jet production is known, and we 
shall summarise its derivation in Section [3.2| . The impact factor for W + 1-jet production 



is derived in Section |3.3| . Next, we summarise jet production at hadron colliders in the 
high-energy limit. 



3 Impact factors 

In order to show how to extract the LO impact factor for W + 1-jet production, we use 
as a paradigm parton-parton scattering, and the derivation of the LO impact factor for 
jet production. 
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yw-(y Jl +y j2 )/2 with | yji -y ja |>2 

Figure 6: Distributions of the rapidity of W boson with respect to the jet average, y v 



(Uji +Z/ja)/2 & t Iz/ji ~Un\ — 2 for the subprocesses of Eq. (|2.2|) at the LHC centre-of-mass 
energy yfs = 14 TeV and with Pj ±m i n = 30 GeV. The dashed line has been generated by 
taking the amplitudes in the high-energy limit, as explained after Eq. (|4.3|) . 



3.1 Dijet production at hadron colliders 

In the high-energy limit s t, the BFKL theory assumes that any scattering process is 
dominated by gluon exchange in the crossed channel, which for a given scattering occurs 
at 0(a 2 s ). This constitutes the leading-order (LO) term of the BFKL resummation. The 
corresponding QCD amplitude factorizes into a gauge-invariant effective amplitude formed 
by two scattering centers, the LO impact factors, connected by the gluon exchanged in 
the crossed channel. The LO impact factors are characteristic of the scattering process 
at hand. The BFKL equation resums then the universal leading-logarithmic (LL) cor- 
rections, of 0{a n s ln n (s/|t|)), to the gluon exchange in the crossed channel. These are 
obtained in the limit of a strong rapidity ordering of the emitted gluon radiation, 

Vl >2/2 > ••• >Z/n-l >2/n- (3.1) 
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For an arbitrary scattering, the LO term of the BFKL resummation is contained in the 
higher-order terms of the expansion in a s . For dijet production in hadron collisions, the 
LO term of the BFKL resummation is already included in the LO term of the expansion 
in a s . In this respect, dijet production in hadron collisions at large rapidity intervals is 
the simplest process to which to apply the BFKL resummation, and thus we shall use it 
as a paradigm. 

Since the cross section for dijet production in the high-energy limit is dominated by 
gluon exchange in the crossed channel, the functional form of the QCD amplitudes for 
gluon-gluon, gluon-quark or quark-quark scattering at LO is the same; they differ only 
by the colour strength in the parton-production vertices. We can then write the cross 
section, 

da = x a f(x a ,/j 2 F )x b f(x b ,fj 2 F )da, (3.2) 
in the following factorised form [E2|, [23|, |2T 



X afeff( x ai Vf) X bfeff( x bi ~nZ~~%Z ' (3-3) 



<Pp a , x <Pp b , x dy a ,dy v aJCJJ " a ™ 0JCJJ " dPp a , x d?p. 
where 

x ° a = , 4 = ^e-y« , (3.4) 



are the parton momentum fractions in the high-energy limit, a' and b' label the forward 
and backward outgoing jet, respectively, and the effective parton distribution functions 
are , 

f eff (x, /4) = G{x, fj.%) + - [Q f (x, /4) + Qf(x, n 2 F )] , (3.5) 

/ 

where the sum is over the quark flavours. In the high-energy limit, the gluon-gluon 
scattering cross section becomes p2| 



dUgg 



C A a s 



d 2 p a ' d 2 p b > Lba'PJ ± ± l\p b . 



f{q ax ,qb ± , Ay) ly-^ , (3.6) 



X 1 



with Ay = y a '—yv and qi x the momenta transferred in the t-channel, with q a± = —p a > ± and 
Qb x = Pv i an< i where Ca = N c = 3. The quantities in square brackets are proportional to 
the impact factors for jet production. We shall analyse them in Section |3.2| . The function 
f{q a± , qb ± , Ay) is the Green's function associated with the gluon exchanged in the crossed 
channel. It is process independent and given in the LL approximation by the solution of 
the BFKL equation. Its analytic form is, 

1 °° f 00 f\a \ 2 \ iv 

f(q a± ,qb ± ,Ay) = ^ c ** / dv e»^* f^U , (3.7) 

- - (27r) 2 |g a± ||g bx | n f-^ J _ x \\q bl \ 2 J 

with the azimuthal angle between q a and q b , and u)(u,n) the eigenvalue of the BFKL 
equation with maximum at u(0,0) = 41n2Cyia s /7r. Thus the solution of the BFKL 
equation resums powers of Ay, and rises with Ay as f(q a± , qb ± , Ay) ~ exp(u;(0, 0)Ay). 
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3.2 Impact factors for jet production 



In the high-energy limit, ln(s/|£|) ^> 1, the LO QCD amplitudes for gluon-gluon, gluon- 
quark and quark-quark scattering all factorise into two impact factors for jet production. 
In order to determine explicitly the impact factors for gluon-jet and for quark-jet produc- 
tion, we shall consider here two of the subprocesses above. 

The cross section for the scattering of two gluon into two gluons, g a Qb 9a' 9b 1 -, at 
LO in the high-energy limit, s = x° a x®s 3> can be written as 



da 



!)!) 



dy a 'd 2 p a > 



An 



■8{Vi4 ~ \Pa>Je y "') 



(3.8) 



x 



dyvd 2 p b > A 

All 



8 2 (Pa> +Pb')\M 



gg^99\ ) 



where \M. 



9 9~*9 9\ 



is the squared tree amplitude, summed (averaged) over final (initial) 
helicities and colours. The amplitude for gluon-gluon scattering, g a g b — ► g a > gy with all 

(3.9) 



external gluons outgoing, can be written as |25|, |26 

'igr' c c"{jp v a %vy 



-\iraabb 
v a v a iv h iv h 



2s 



igf Wc C"{f b \pl) 



(3.10) 



where the z/s label the helicities and the LO vertices g* g — > g are given by 

P\ 
Pb' ± ' 

where we represent the transverse momentum p±_ on the complex plane, p± = p x +ip y . The 
functions C transform into their complex conjugates under helicity reversal, C*({k u }) = 
C({k~ u }). The helicity-flip function C(p + ,p f+ ) is subleading in the high-energy limit. 
Thus each function, C", has only the two helicity configurations of Eq. ( |3.10| ) allowed in 
the high-energy limit. 

We define the impact factor, I(p a ,Pa') or I(Pb,Pv), as the square of each term in 
squared brackets in Eq. ( |3.9|) , summed (averaged) over final (initial) helicities and colours. 
Thus the LO impact factor for a gluon jet is 

1 in /"•'•' ' ( " ; ''f ,, 

3 ir r 



I 9 (Pj,Pj>) 



2(iV c 2 - 1) 



JV*-1 

with j = a,b and with implicit sums over repeated indices. 

From Eq. ( ft.9|) and ( |3.11| ), the squared amplitude for gluon-gluon scattering is 

|2 



(3.11) 



\M 



99^991 



As 2 

-J^I 9 (Pa,Pa>)I 9 (Pb,Pb> 

AC 2 



N*-l y t 2 



9 4 s 2 
- q — 
2 y t 2 



(3.12) 
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Analogously, the quark-gluon q a — ► q a i gy scattering amplitude in the high-energy limit 
is 



njaa'bb' 
m v a v h v h i 



2s fcA^C*(p* jtf*)] j Ugf bb ' c C°\ Pb »,pir 



t L 



(3.13) 



with LO vertices g* q — > q, 



1/2 



(3.14) 



Again, each function, C qq , has two helicity configurations allowed. The LO impact factor 
for a quark jet isQ 



I q (Pa,Pa> 



2N r 



2N r 



[gK'- a c qq {p v a,p7 a )\ 9>4a> [c qq {p v a \P7 a )T 



(3.15) 



The squared amplitude for quark-gluon scattering is then 



\M I 2 



4s 2 

-£rI q (Pa,Pa>)I 9 (j)b,Pbi 



2C 



A „4 



Nr. 



g-- = 2g^ 



t 2 



t 2 



(3.16) 



Since Eqs.( p.l2"|) and (|3.16|) only differ by the colour factor 9/4, to calculate the dijet pro- 
duction rate, Eq. fl3.3|) , it suffices to consider one of them, such as gluon-gluon scattering, 
and include the others through the effective p.d.f. in Eq. (|3.5| ). Using Eq. ( |3.12|) and 
replacing t 2 — > \p 2 a , ||p 2 , |, the cross section for gluon-gluon scattering (|3.8|) becomes 



d 2 p a > d 2 p b > 



\Pa' I' 



-$ 2 (Pa\ +Pb\ 



\Pb\ I 2 



(3.17) 



At higher orders, powers of ln(s/|t|) arise, which can be resummed to all orders in 
a s ln(s/|t|) through the BFKL equation. The term between square brackets in Eq. ( [3.171) 
is the LO term of the BFKL resummation. Since in performing the resummation the 
factorization formula ( |3.3D holds unchanged, to obtain the resummed parton cross section 
it suffices to replace the LO term of the BFKL ladder in Eq. (|3.17|) with the full ladder, 
thus obtaining Eq. (|3.6|) . 



3.3 The impact factor for W + 1-jet production 

In W + 2-jet production in the limit [y^ — yj 2 \ ^> 1 (see Section |2.2|) , the parton subpro- 
cesses qq — > Wqq and qg — > Wqg and qq — > Wqq all feature gluon exchange in the 

^We use the standard normalization of the SU(iV c ) matrices, tr(A c A c ) = S cc /2. 
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crossed channel. Thus the functional form of the corresponding QCD amplitudes is the 
same. They differ only by the colour strength in the impact factor for jet production, sep- 
arated from the impact factor for W boson and the other jet by the gluon in the crossed 
channel. It suffices then to consider only one of the subprocesses above. We shall take 
qg — > W qg and we shall suppose, for the sake of clarity, that the W boson is produced 
in the positive-rapidity hemisphere. The subprocesses above factorise according to the 
kinematics 



Vv 



Vq > y g 



\Pw ± I - \p 



\Pg 



In the high-energy limit, the cross section for qg — > W qg scattering is 



da 



wqg 



dy w d 2 Pw dy q d 2 Pq± , _ _ , 

47r(2vr)2 4vr °\V 8X ° m±e > 



(3.18) 



(3.19) 



$ 2 (Pw ± + Pq ± + PgJ \M qg ^ W qg\ 2 , 



with the W transverse mass as in Eq. (|2.3|) . If we include the subsequent decay of the W 
boson into a lepton pair, the kinematics in the high-energy limit becomes 



Vq > Vg 



\Pe 



\Pu ± I - \P 



\P. 



9x1 



(3.20) 



The cross section for qg — > qg(W — >)eu scattering can be obtained from Eq. ( |3.19| ). 
replacing the W boson with the lepton pair and using the amplitudes calculated in 
Refs. [^7|, . In the notation of Ref . 0] the colour decomposed amplitude is 

A(l g ,2,3,4 5 ;5 e -,6 e ) (3.21) 

= gW^wism) ^(A^)A a ^))j4 A 6 (l q ,a(2),a(3),4 q ;5 B ,6 e ), 
aes 2 

where legs 1, 4 are the qq pair, legs 2, 3 are the gluon legs, and legs 5,6 are the lepton pair; 
gw is the weak coupling and Vw(s) the W propagator 



V w (s) 



s-M^ + iTwMw 



(3.22) 



where Mw and Tw are the mass and width of the W. The colour ordered subamplitudes 
are 

2 



A 6 (l-,2-,3 + ,4+5t,6,- 



S23 L 



[13] (21) [5 4] (6|(1 + 2)|3> 



[21]* 



123 



[4 3] (2 4) (6 1) (2 1 (3 + 4) 1 5) (2 1 (3 + 4) 1 5) (6 1 (1 + 2) 1 3) 



(4 3)t 23 4 



A 6 (l-,2 + ,3-,4+,5±,67) = — 



2 

^23 L 



[2 1] (4 3) 
(31) 2 [54] (6|(1 + 3)|2) 



+ 



(21)* 12 3 

[42] 2 (61)(3|(2 + 4)|5) , (3 1) [4 2] (6 1) [5 4] 



[4 3] t 2 34 



(21) [4 3] 



(3.23) 
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with the spinor products defined in Appendix [A]. The subamplitudes ( |3.23|) are symmetric 
under the exchange [|30] 



4, 2 



6 , (ij) <-> [ji] ■ 



(3.24) 



In Eq. ( p.23|) we have neglected the subamplitudes with like-helicity gluons because they 
are subleading in the high-energy limit. 

Next, we make the correspondence p^ = p a , p\ = p a >, P2 = Pb and p% = py and 
according to Eq. ( |3.23| ) we always identify (5)6 as the (anti)lepton momentum. The 
amplitude for q a g b — > q a i g v {W — >)ev scattering ( |3.21| ) in the high-energy limit, y q ~ 
y e ~ y u ^ VVi is obtained by computing the sub-amplitudes ( |3.23|) in the corresponding 
kinematics (Appendix ||) 



A 



ia gb->q a ' 9b' e v 



(3.25) 



2S [gK>- a C« q (f a \P7\Pe,Vu,q)\ j [igf Wc C^{jpl\plY 



with py = Pg, and q± the momentum transferred in the crossed channel, q± = py ± , and 
t ~ — | | 2 , and the vertex g* d — > u e~ v given by 



C du {p+,p q ,p e ,Pu,q) 



igwVw(Sve) 

yj {Pq + Pw)pt 



(3.26) 



x 



(PqPe) 



PUl - (Pq + Pw) Ptx iPt * Pw±Pt + Pt Pe ± 



tabb' 



pi 



ta'bb' 



with 



ta'W 
tabb' 



(Pa +Pw) 2 ^ ~Pq P W ~ \Pw± I 
(Pq+Pw) 2 , 



(3.27) 



and with p w = p e + p v - In the argument of the vertex in the left-hand side of Eq. ( |3.26| ) 
we do not write explicitly the helicity of the lepton pair, since that is uniquely fixed by 
the helicity of the quark pair and of the W boson. The vertex g* u — > d e + v is obtained 
from Eq. ( |3.26|) by exchanging (p„ <-> p e ) 

C ad (p+,p-,p e ,Pv,q) = C* u (p+,p q ,p u ,p e ,q). (3.28) 

Using the symmetry ( ft.24 ) of the vertex (|3.26| ), we obtain the vertex g* d — > ue + v 

C da (p-,p+, Pe ,Pv,q) = - [c du (pi,p-,p e ,Pu,q)]* ■ (3.29) 
The vertex g* u — > d e~ v is then obtained from Eq. ( |3.29| ) by exchanging (p u «-» p e ) 

C u \p-,pi,Pe,Pv,q) = C du (p-,p+ Pl/ ,p e ,q). (3.30) 
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The impact factor for W + 1-jet production, I qw , can be obtained by squaring any of 
the effective vertices (|3.26| - |3~30D and by integrating out the lepton pair; however, by using 
Eqs. ( |3.21| ) and (|3.23| ) we have computed directly the squared amplitude for q g — > q g W 
scattering, and compared it to Ref. Taking then the high-energy limit ( 3.18|) , the 

squared amplitude summed (averaged) over final (initial) colours and helicities, reduces 
to 



\M 



qg->wqg\ 



4S 2 



I qW (Pa,Pq,Pw,q) I 9 (j>b , Pb> ) 



with 

I 9W (Pa,P q ,Pw,q) 



Xcc' 2 

g 2 

ZNctabb'ta'bb' 2 



m. 



tabb' 



Z + 



tg'W 
tabb' 



i (i + z 2 ) 



where we have defined the momentum fraction 



z = 



pi 



p+ + Pv 



(3.31) 



(3.32) 



(3.33) 



Using Eq. ( |3.27| ) and i ~ — \q±\ 2 , we can rewrite the impact factor ( |3.32| ) as, 



I qw (Pa,P q ,Pw,q) 

6 cc' 



29 w 



m 



2 ~z\q±\ 2 + \p q± \ 2 - \P wA 2 f + |^| 2 , x + ^ 



w 



tabb'ta'W 



(3.34) 



In the small \q±\ limit, the jet opposite to the impact factor for W + 1-jet production 
becomes collinear, and the cross section obtained from the squared amplitude ( |3.31| ) yields 
an infrared real correction. Since the latter may have at most a logarithmic enhancement 
as — > 0, the squared amplitude (|3.31|) cannot diverge more rapidly than l/|g^| 2 . This 
entails that in the small \q±\ limit, the impact factor ( |3.34|) must be at least quadratic in 
|g_ij, I qw ~ (9(|gjJ 2 ). Using q± = —(p q± +p w± ), it is immediate to see that this is the 
case for Eq. (|3.34j) . In addition, as q± — > we have an almost on-shell gluon scattering 
with a quark, then p q± — > —p w± and averaging over the azimuthal angle of q±, Eq. ( ft.34j ) 
becomes 



lim I qw 



S cc ' a 2 
2N r 2 



\q±\ z(l - z) 
\pwA 2 + zm w 



{l + z 2 )(\p w± 


4 + z 2 ™^) + 4z 2 m 2 v 


Pw ± 


2 


z{\Pw ± 


2 + zm\) 2 



45 c 



\q±\ z(l - z) 



| 2 + zm% 



\M 



qg^wq\ j 



(3.35) 



which explicitly shows that the impact factor is positive definite and that it factorises into 
the squared amplitude for q g — > W q scattering [^(J , as it should. 
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4 The production rate for W + 2 jets 



In the collision of two hadrons A and B, the differential production rate of a If boson 
with two associated jets is given in terms of the rapidities and transverse momenta by 



da 



d 2 p jl± d 2 p h± d 2 p w± dy h dy h dy v 



X afi/A(Xa, /4) %bfj/B(Xb, /4) ^ ( Pw ^ + Ph± + Pj 



(4.1) 



with parton momentum fractions ( |2.3|) . In Eq. 
fully contained in the squared amplitude. 



J]) the dynamics of the scattering are 



In the limit \y^ — yj 2 \ ^> 1, as discussed in Section |3l| we can identify an outgoing 
parton with a (anti) quark, while the other, that we called b' according to the notation of 



Section |3]3|, can either be a quark or a gluon. The cross section for W + 2-jet production 
can be written in the factorised form (|3.3|), by substituting Eq. ( |3.31| ) and using Eq. 



da 



d 2 p q± d 2 p b > ± d 2 p w± dy q dy v dy w 



jqwjg 



Qbj. 



327r 5 |g a J 2 |g fc± | 2 



(4.2) 



where g a± = — p q± —p Wl _ and q b± = Pb' ± , and where we have substituted t 2 with \q a± \ 2 \qb ± \ 2 - 
In the first p.d.f. the sum is over (anti)quark flavours, and the impact factors are given in 
Eqs. ( 3.11 ) and (|3.32 ). The last term is the LO term of the BFKL resummation. Thus, 
to obtain the BFKL-resummed cross section we just need to replace S 2 (q a± — q b± )/2 with 
f(q a± ,qb ± ,Ay), as in Eq. ( ETTp - 

However, in Eq. (|4.2[) energy and longitudinal momentum are not conserved. The 
parton momentum fractions in the high-energy limit, x° a and x®, given in Eq. (|3.19| ) un- 
derestimate the exact ones ( j2.3| ) and accordingly the p.d.f. 's can be grossly overestimated. 
Thus for numerical applications and for a comparison with experimental data, it can be 
convenient to perform the high-energy limit only on the dynamical part of Eq. ( |4.l| ), by 
writing the squared amplitude in the factorised form (|3.31|) , while leaving the kinematics 
untouched. This entails 



da 



d 2 p q± d 2 p b > d 2 p w± dy q dy b/ dy v 



2^ %aQi(%a, /J>f) x bfeff(%b, ^f)'^' 



32vr £ 



jqwjg 
t 2 



S 2 (q a± - qb s 



(4.3) 



For the invariants t and t a ' b b', implicit in the square brackets, two options are possible:| 
■^The invariant t a bb> is the same in the exact and high-energy kinematics. 
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Figure 7: The VF + 2-jet production rate as a function of the rapidity interval between the 
jets Ay = \yj 1 — y 3 - a |. The solid curve is the exact production rate (|4.1| ); the dot-dashed 
curve is the production rate in the high-energy limit (O); the two dashed curves are given 
by the production rate (|4.3|) , with the two approximations for the dynamics mentioned in 
the text. 



(a) they are taken to be exact, namely t = 2pb-py and t a i h y = (p a i +p w ) 2 . For instance, 
the dashed lines of Fig. ^| have been obtained from Eq. Q4.3| ) with option (a); 



(b) t and t a >bb' are in the high-energy limit, as defined below Eq. fl3.25p and in Eq. ( |3.27| ) 
respectively. 



Note that Eq. (|4.3|) , with the two approximations for the dynamics above, and Eq. (|4.2j ) 
have the same theoretical validity, however their numerics may be rather different. In 
order to examine that in detail, in Fig. [F] we consider W + 2-jet production as a function 
of the rapidity interval between the jets Ay = \yj 1 — yj 2 \. For the renormalisation and 
the factorisation scales we keep the same choice as in Section |2.2| . The solid curve is 
the exact production rate ( |4.1| ); the dot-dashed curve is the production rate in the high- 
energy limit ( |4.2| ); the two dashed curves are given by the production rate ( |4.3|) , with 
the two approximations listed above: (b) is the upper dashed curve, and (a) is the lower 
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one. Note that the exact production rate is contained between curves (a) and (b), with 
(b) yielding the best numerical approximation to the exact curve, while the high-energy 
limit on the LO kinematics and dynamics (the dot-dashed curve) is rather distant from 
the exact production rate, unless Ay is quite large. The range between curves (a) and (b) 
may be viewed as a band of uncertainty on the high-energy limit at LO. 



5 The BFKL Monte Carlo 

In a comparison with experimental data, it must be remembered that the LL BFKL resum- 
mation makes some approximations which, even though formally subleading, can be nu- 
merically important: a) The BFKL resummation is performed at fixed coupling constant, 
thus any variation in the scale at which a s is evaluated appears in the next-to-leading- 
logarithmic (NLL) terms, b) Because of the strong rapidity ordering any two-parton 
invariant mass is large. Thus there are no collinear divergences in the LL resummation 
in the BFKL ladder; jets are determined only at tree-level and accordingly have no non- 
trivial structure, c) Finally, energy and longitudinal momentum are not conserved, since 
the momentum fraction x of the incoming parton is reconstructed from the kinematic 
variables of the outgoing partons only, and not including the radiation from the BFKL 
ladder. Therefore, the BFKL theory can severely underestimate the exact value of the 
x's, and thus grossly overestimate the parton luminosities. In fact, if a W boson + (n + 2) 
partons are produced, we have 

Xa(b) = \PllA e nvn + \P^A e (-)y> 2 + !^ e (-W + V *£ e (-)w (5.1) 
V s V s V s j=iv S 

where the minus sign in the exponentials of the right-hand side applies to the subscript b 
on the left-hand side. In the BFKL theory, the LL approximation and the kinematics (|3.1| ) 
imply that in the determination of x a (x^) only the first (last) term in Eq. (|5.1| ) is kept. 
The terms neglected in Eq.( |3.4| ) are formally subleading. However, a comparison of three- 
parton production with the exact kinematics ( |5.1|) to the truncation of the BFKL ladder 
to 0(a^) shows that the LL approximation entails sizable violations of energy-momentum 
conservation p2| . 

Kinematic cuts and constraints like Eq. (|5.1f) can be implemented in the BFKL frame- 
work by unfolding the BFKL integral equation, resulting in an explicit sum over the 
number of emitted gluons. Each term in this sum is then a phase space integral over the 



BFKL gluon phase space, which allows for a BFKL Monte Carlo to be constructed [33, 34]. 
Since each emitted BFKL gluon enters the calculation with an explicit phase space inte- 
gral, this approach allows for the running of the coupling to be included into the BFKL 
solution, and also for the gluon radiation to be taken into account in Eq. (|5.1|). 



The first step in this procedure is to transform the relevant Green's function f(q a ±, qb±, Ay) 
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of Eq. (p.7|) to moment space via 

f(q a ±,qb±,Ay) = J ^e wA ^(<f a± ,<f 6± ) . 
We can then write the BFKL equation as 



L 7T 

where the kernel K is given by 

— * — * 

K(q a ±, q b ±, k±) = f u [<fa± + k±, q b± ) 



a* f d 2 k i 



K(q a± ,q b± , kj 



k\ + (q a ± + k A 



fui{q a ±, q b ±] 



(5.2) 



(5.3) 



(5.4) 



and a s = ol s -^. The first term in the kernel accounts for the emission of a real gluon of 
transverse momentum k± and the second term accounts for the virtual radiative correc- 
tions. 

We now separate the k± integral in ( |5.$ ) into 'resolved' and 'unresolved' contributions, 
according to whether \k±\ lie above or below a small transverse energy scale fi. The scale 
/i is assumed to be small compared to the other relevant scales in the problem (such as the 
minimum transverse momentum, for instance). The virtual and unresolved contributions 
are then combined into a single, finite integral. The BFKL equation becomes 



w f u (q a ± ,ql±) = -5 (q a± - q b± 



d 2 k\ 



n Jk^>u? k 



+ — 



a s [ d 2 k 



7T / k , 



2 U(q a ± + k±,q b ± 
fu>(q a x + k±, q b ±) 9(n 2 - k]_) 



ql± fu(qa±i gbjj 
k\ + {q a ± + k ± ) 



(5.5) 



The combined unresolved/ virtual integral can be simplified by noting that since k\ <C 
1a±y c&i. by construction, the k± term in the argument of f w can be neglected, giving 



1 a s 
(v - uo) fu(q a ±, qb±) = ^(ii- Qb±) H S ~ 



d 2 k A 
~kj 



where 



a. f d 2 k i 



u> Q = — 

7T 



e^ 2 - k\) 



kl + (q a± + k ± ) 2 



— * 

fu{&± + k±,q b ±) , (5.6) 



a s In ( ) . (5.7) 



The virtual and unresolved contributions are now contained in u Q and we are left with an 
integral over resolved real gluons. Eq. ( |5.6|) is now solved iteratively, and performing the 
inverse transform we have 



f(q a ±,q b ±,&y) = ^2f {n) (q a ±,qb±, Ay) 



(5. 



n=0 



19 



where 



f {0 \q a ±,q b x,Ay) 

f {n - 1] (q a ±A b ±^y) 



U 2 ~ 


a s Ay 


k\. 




U 2 ~ 


a s Ay 


a 2 
. ( ia±. 









Qb±, 



Y[ / d 2 &i± % ^i? - S(q a± -q b ±-^2 ki±) , 

1=1 ^ J 8=1 

(g a ± + J2Z\ kj±) 2 



e(kf ± - o(y l -i - % 



{q a ± + ELi k j±) 



. (5.9) 



Thus the solution to the BFKL equation is recast in terms of phase space integrals for 
resolved gluon emissions, with form factors representing the net effect of unresolved and 
virtual emissions. In this way, each depends on the resolution parameter /i, whereas 
the full sum / does not. 

The derivation given above only applies for fixed coupling because we have left a s 
outside the integrals. The modifications necessary to account for a running coupling 
a s {kf L ) are straightforward [ 33| . In the rest of this paper, we will however discuss only 
the fixed coupling version of the BFKL Monte Carlo with the coupling entering the BFKL 
equation set to a s (p 2 ±m i n ), and with energy momentum conservation built in through 
Eq. ( p.l| ). The effects of including the BFKL gluon radiation in the Bjorken x's are far 
bigger than the effects of the running coupling, which amount to an approximately 10% 
effect in the cross section of Fig. H. 



6 BFKL observables 

After having considered several approximations to the high-energy limit and introduced 
the BFKL Monte Carlo as the tool that we shall use to analyse the BFKL gluon radiation, 
we turn now to the analysis of the effects of the BFKL radiation on some physical observ- 



ables. From Eq. (|3.7|) we see that in order to detect evidence of a BFKL-type behaviour 
in a scattering process, we need to obtain Ay as large as possible. In the context of dijet 
production this can be done by minimizing the jet transverse energy, and maximizing the 
parton centre-of-mass energy s. Since s = x a Xf,s, in a fixed-energy collider this is achieved 
by increasing the parton momentum fractions x at b, and then measuring the dijet produc- 
tion rate da/dAy. However, in dijet production three effects hinder the characteristic 
growth of the BFKL ladder ( |3.7D with respect to LO production: 



as the x's grow the parton luminosities fall off, making it difficult to disentangle the 



eventual BFKL-driven rise of the parton cross section from the p.d.f.'s fall off |f23|. Q 



the implementation of the exact x's (|5.1| ) in the BFKL Monte Carlo |35|], rather 



than using x° a b ( |3.4| ) as prescribed by the high-energy limit, shifts the p.d.f.'s toward 
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Figure 8: The W + 2-jet production rate as a function of the rapidity interval between 
the jets Ay = \y^ — yj 2 \, with acceptance cuts y w , yj 2 > 1 and y^ < — 1, or y w , yj 2 < — 1 
and > 1. The diamonds are the exact production rate fl4.1|) ; the dashed curve is the 
production rate in the high-energy limit ( |4.3|) with option (a); the dotted curve is the 
same with option (b); the solid curve includes the BFKL corrections. 



smaller values, and thus further suppresses the production rate. This effect is already 
present at 0(a\) H; 



in dijet production both the tagged jets have typically the same minimum transverse 
energy; at NLO, the dijet cross section as a function of the difference T> between 
the minimum transverse energies of the two jets turns out to have a slope da/dD 
which is infinite at T> = p6|, B7I]. This hints to the presence of large logarithms 



of Sudakov type, which can conceal the logarithms of type ln(s/t) characteristic of 
the BFKL dynamics.^ 



^Logarithms of Sudakov type are contained in the BFKL solution (3.7), however they lack the running 
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The combination of these three effects changes drastically the shape of dijet production 
da/dAy as a function of the rapidity interval between the tagged jets, showing a deple- 
tion rather than the characteristic increase of the BFKL analytic solution. 

In Fig. | we consider W + 2-jet production as a function of Ay, and with acceptance 
cuts y w , yj 2 > 1 and y^ < — 1, or y w , yj 2 < —1 and > 1. For all of the curves of 
Figs. |8|-|I2| we choose fi R1 = pj 1± and \i m = (pj 2± + m±)/2 as renormalisation scales, and 



/i Fi = /i F2 = (biixl + \Ph± \ + m J-)/2 as factorisation scales. We justify the peculiar scale 
choices above as follows: we note that our calculations are at LO (from the renormalisation 
point of view), thus the scale choice is completely arbitrary, as long as it is physically 
unambiguous. However, a uniform choice for all of the curves in the same figure is required 
for a consistent comparison between different approximations. In addition, in the high- 
energy limit the impact factors for iy + 1-jet production on one side and for jet production 
on the other can be viewed as two almost independent scattering centres linked by a gluon 
exchanged in the crossed channel, thus it makes sense to run a s according to the scale 
set by each impact factor. Accordingly, in the LO calculation ot 2 s must be understood 
as a s (p 2 i± ) a s ((pj 2± + rn±) 2 /4:). In the high-energy limit it is possible (and would make 
sense) to choose the factorisation scales equal to the renormalisation scales, however for 
the exact production rate this choice would not be physically sensible since no high-energy 
factorisation is present, thus for the factorisation scales we keep the same choice as in the 
previous figures. In Fig. |] the diamonds represent the exact production rate (f4.1|) ; the 
dashed curve is the production rate in the high-energy limit ( |4.3| ) with option (a); the 
dotted curve is the same with option (b); the solid curve includes the BFKL corrections. 
In Figs. p|-|T2"| we have computed the BFKL corrections using Eq. ( f4.3| ) with option (b). 
However, the particular option we choose is immaterial since the uncertainty related to 
the choice of option in Eq. ( |4.3p is much smaller than the uncertainties intrinsic to the 
BFKL resummation, the latter being due to the leading-log approximation, the choice of 
scale of a s and the approximation on the incoming parton momentum fractions. Note 
that the curve of Fig. |8] is both qualitatively and quantitatively different from da/dAy 
in dijet production: the peak in Fig. 13 is a striking confirmation of the dominance of the 



configurations asymmetric in rapidity, discussed in Section 2.2. In fact the symmetric 



acceptance cut strongly penalises the asymmetric configurations when Ay approaches its 
minimum value; since the asymmetric configurations dominate the W + 2-jet production 
rate, the effect is a strong depletion of the latter. In addition, the BFKL ladder (solid 
curve), which includes energy-momentum conservation ( |5.1|) , shows a substantial increase 
of the cross section with respect to the LO analysis (dotted and dashed curves), as opposed 
to a decrease in the dijet case. To understand how this comes about, we note that the 
presence of at least three particles in the final state makes the threshold configurations, 
and thus the logarithms of Sudakov type, much less compelling than in the dijet case. 
Secondly, the implementation of the kinematic constraint (|5.1|) in the BFKL Monte Carlo, 



of a s and they are not consistently resummed. 
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Figure 9: The ratio f(x°, /ip)/f(x, of the p.d.f. as a function of the x°'s ( p.!9| j 
in the high-energy limit versus the p.d.f. as a function of the exact x (|5.1|) ; the 
dashed-dotted curve is the ratio f{x G a , /i' F ) / f(x a , /ip) , and the solid curve is the ratio 

feff(x° b , l4)/ feff(%b, Hi)- 



rather than using x Q ab ( |3.19| ) in the high-energy limit, has a much lesser impact than in 
the dijet case. This is due to the fact that the valence quark distribution in q g — > q g W is 
much less sensitive to x variations than the gluon distribution in g g — *■ g g. To analyse this 
more precisely, we consider in Fig. || the ratio f(x°, up) / f(x, fip) of the p.d.f. as a function 
of the x°'s (|3.19|) in the high-energy limit versus the p.d.f. as a function of the exact x (|5.1| ). 
The ratio is calculated for each event in the Monte Carlo as the ratio of the p.d.f. evaluated 
at x° compared to an evaluation at x, weighted with the contribution of this event to the 
cross section according to (|4.3|) with option (6) and the BFKL ladder added. Finally, this 
distribution is binned in Ay. To be definite, since the high-energy factorisation entails that 
each impact factor is associated to one of the two incoming partons, we can term the ratio 
f(x®, ftp)/ f(x a , }J, F ) as the one associated to the impact factor for W + 1-jet production, 
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and the ratio f e ff(xl,fj, 2 F )/f e ff(xb,f^) as the one associated to the impact factor for jet 
production. As we see from Fig. ^], the solid curve is much farther away from 1 than the 
dashed-dotted curve. Since the effective p.d.f. is dominated by the gluon distribution, 
this implies that the ratio f e ff(%bi A*f)/ feff( x b, /40 * s m uch more sensitive to variations of 
the x's than the ratio /i|)//(x a , //§,), which is made by valence quark distributions. 
Accordingly, we obtain a smaller depletion of the BFKL Monte Carlo prediction in W + 2- 
jet production as compared to dijet production. In addition, both the curves in Fig. |9| rise 
as Ay grows. That entails that the BFKL radiation, which enters the determination of 
the x's in the denominator, yields as expected a contribution which is growing with Ay. 




Figure 10: The average azimuthal angle (cos A<p), where A0 



'31 



'32 



— 7r, as a function 



of the rapidity interval between the jets Ay = \yj 1 —yj 2 \, with acceptance cuts y w , yj 2 > 1 
and yj 1 < —1, oi y W) yj 2 < —1 and y^ > 1. The diamonds are the exact production rate 
( O) ; the dashed-dotted curve is the production rate in the high-energy limit fl4.3Q with 



option (6); the solid curve includes the BFKL corrections. 
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A variable that has been extensively studied as possibly sensitive to BFKL effects is 
the azimuthal angle decorrelation A0 = \<pj 1 — <pj 2 \ — it between the most forward and 
backward jets in inclusive dijet samples. At LO the jets are supposed to be back to 
back, with a correlation which is smeared by gluon radiation induced by parton show- 
ers and hadronization. However, if we look at the correlation also as a function of Ay, 
we expect the gluon radiation between the jets to further blur the information on the 
mutual position in transverse momentum space, and thus the decorrelation to grow with 
Ay. From Eq. (|3.7|) , we see that the BFKL-induced gluon radiation might account for 



that |23|, |4], |32], |33|, |38|| . The decorrelation between the tagging jets has been analysed, 
and indeed observed, by the DO Collaboration in dijet production at the Tevatron Col- 
lider However, the BFKL-induced radiation predicts a stronger decorrelation than 
the data, even though the BFKL Monte Carlo |33| shows a much more realistic azimuthal 
decorrelation than the BFKL analytic solution. The data are correctly reproduced by 
the HERWIG Monte Carlo generator [IM HO, ET|, which includes parton showers and 



hadronization. This suggests that the azimuthal angle decorrelation Acj), picking up pref- 
erentially configurations where the tagged jets are back to back, is sensitive to threshold 
configurations, and thus to logarithms of Sudakov type, even more than it is in the inclu- 
sive dijet production rate da/dAy |37|, However, as discussed in the paragraph above, 
in W + 2-jet production we expect the logarithms of Sudakov type to play a much less 
significant role. Thus, in analogy with dijet production, in Fig. 10 we consider the average 
azimuthal angle (cos A<p) as a function of the rapidity interval between the jets Ay. The 
acceptance cuts are the same as for Fig. H The diamonds are the exact production rate 



(4.1); the dashed-dotted curve is the production rate in the high-energy limit fl4.3|) with 
option (b); the solid curve includes the BFKL corrections. The average azimuthal angle 
being defined as a ratio of production rates is much less sensitive to scale variations than 
the curves of Fig. §. 



In Fig. [TT| we plot the BFKL prediction for the mean number of jets (n) with p± > 
Pj±mm = 30 GeV, emitted by the BFKL ladder as a function of the rapidity interval 
between the jets Ay, and the BFKL prediction for the same variable in the rapidity range 
— 1 < y 3 ; < 1. We see that the mean number of jets rises approximately linearly with Ay, 
and accordingly that the mean number of jets in the rapidity range — 1 < yj < 1 stays 
constant. We can crudely understand this, by noting that for a very large Ay the cross 
section from Eqs. Q3.7| ) and ( |4.2j ) behaves like 

(w(0,0)Ay) n 



i „uj(0,0)Ay 

" ni 

n=0 



^ > ^ 



with u(0,0) = 41n2Cyia s /7r, and a power of a s for each real correction to, and therefore 
for each emitted gluon from, the BFKL ladder. Up to corrections of type ln(p_|_/p,,_|_ min ) 



34] , the mean number of jets emitted by the BFKL ladder is then 



(n) = [na . )lAy ~Lu(0,0)Ay. (6.2) 
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Figure 11: Solid: the BFKL prediction for the mean number of jets (n) withpx > Pj±mm = 
30 GeV as a function of the rapidity interval between the jets Ay — \yj\ — Uj2 

I. Dotted: 

the same in the rapidity range — 1 < yj < 1. 



For a s (p|j_ min ) with Pj± m i Il = 30 GeV, this yields typically a jet each second unit of 
rapidity, which is in rough agreement with Fig. |LT[ 

Finally, in Fig. ^ we consider W + 2-jet production as a function of the transverse 
momentum q a± = — (p q± + Pw ± ) exiting from the impact factor I qw for W + 1-jet produc- 
tion. At LO, q a± = qb ± = Pb' ± , thus q a± is bound to be equal to the transverse momentum 
of the jet opposite to the impact factor I qw (and thus to be always larger than 30 GeV). 
In presence of the gluon radiation of the BFKL ladder, this is not longer true, and q a± is 
allowed to go to zero. However, a simple power-counting argument shows that the pro- 
duction rate is finite as q a± — ► 0. In fact from Eq. ( |3.34|) we know that I qw ~ 0(\q a± \ 2 ). 
Substituting it, and the ladder ( |3.7| ) which goes like 0(l/\q a± \), in Eq. ( |4.2|) , we see that 
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Figure 12: The W + 2-jet production rate, including BFKL corrections, as a function of 
the transverse momentum q a± , with q a± = — (p q± + p w± )- 



as far as the behaviour in q a± is concerned, 



7 Conclusions 

In Section || we have examined the exact LO inclusive rapidity distribution for W + 1- 
jet and W + 2-jet production; as for the latter, we have seen that the dominant parton 
sub-process qg — > qgW produces a great deal of W bosons forward in rapidity. This is 
due to the different shape of the p.d.f.'s of the incoming quark and gluon, and to gluon 



f^ 2 (?a ± +Pq ± +Pw ± ), 



(6.3) 



and therefore the distribution da/dq a± is finite as q, 



0, in agreement with Fig. [12. 
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exchange in the crossed channel which loosens the bound between the W boson and a jet 
on one (rapidity) side, and the other jet on the other side. In Section ^ we have derived the 
impact factor for W + l-jet production, both as a function of the P^-boson momentum and 
including the leptonic decay of the W boson. In Section |] we have compared several high- 
energy approximations at LO to the exact production rate. The range between the most 
extreme high-energy approximations may be considered as the theoretical uncertainty on 
the high-energy limit at LO. 

In Section |6] we have considered some BFKL footprints, most notably the rate da/dAy 
and the azimuthal angle decorrelation da/dAcf) as functions of the rapidity interval Ay 
between two tagged jets. These observables had already been considered in inclusive 
dijet production, however because of the dominance of the configurations asymmetric 
in rapidity and the presence of at least three particles in the final state, which makes 
threshold configurations less relevant, in W + 2-jet production da/dAy and da/A(p take 
on a completely new light. In addition, we have considered the mean number of jets, 
which as expected rises approximately linearly with Ay. Finally, we have computed the 
transverse momentum distribution of the impact factor for W + 1-jet production. At LO 
this is bound from below by momentum conservation at the minimum transverse energy 
of the jet opposite to the W + 1-jet configuration, but with additional gluon radiation it 
is allowed to reach zero, where the distribution is finite. 

Finally, we note that one of the leading contributions to the WW + 2-jet production 
rate in the high-energy limit is obtained by convoluting two impact factors for W + 1-jet 
production with a gluon exchanged in the crossed channel. The analysis of this process, 
in the exact case, in the high-energy limit and with BFKL corrections is left for future 
investigations |H3 . 
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A Multiparton kinematics 

We consider the production of a W(— > ev) and two jets of momenta p a > and py, in the 
scattering between two partons of momenta p a and pb (p° a < and p° < 0).Q 

^Conventionally, in the hclicity amplitudes all momenta are always taken as outgoing. Partons which 
in a physical channel are incoming are then identified by the (negative) sign of their energy. 
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Using light-cone coordinates p = p ±p z , and complex transverse coordinates p± = 
px _|_ jpy^ w ith sca l ar product 2p ■ q = p + q~ + p~q + — p±q]_ — p* ± q±, the 4-momenta are, 

Pa = (p+/2,0,0,p+/2) = (p+ 0;0,0) , 

Pb = (p 6 -/2,0,0,-p 6 72) = (0,p 6 -;0,0), (A.l) 

V l = {(pi +Pi)/2,Re[p l± ],lm[p t± ],(pt -pr)/2) 

= (\Pi ± \e m ,\p i± \e~ m ;\p i± \ cos (f)i,\p i± \ sin <f)i) , % = 3(6'), 4(a), 5(f), 6(e) 

where the first notation is the standard representation p M = (p°,p x ,p y ,p z ), while in the 
second we have the + and - components to the left of the semicolon , and to the right 
the transverse components, y is the parton rapidity and is the azimuthal angle between 
the vector p± and an arbitrary vector in the transverse plane. From the momentum 
conservation (i = b', a', 1, £), 

6 

= ^Pi ± , 

6 

Pa = ~Y.Pt, (A.2) 

6 

Pb = ~J2 P i' 

i=3 

the Mandelstam invariants may be written as, 

= 2pi ■ Pj = pfpj + Pipj - p i± P* j± - P* ± p j± ■ 

so that 

6 

s = 2p a -p b = Y^ pfpj , 

6 

Sai = tya ■ Pi = -^PiPj ' ( A - 3 ) 

i=3 

6 

s M = 2p b ■ pi = - Pfpj ■ 

3=3 



The spinor products are defined as 



(Pi - \Pj+) = (u) , 

(Pi + bi-) = M , (A.4) 

(Pi - \Mpj-) = ■ 
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Using the above spinor representation, the spinor products for the momenta (|A.l ) are 




Pj Pi 
[PiPj) = Pi ± \l ~ Pj ± \l —+ , 

Pi 




-T) 

(PaPi) = —Pixi ( A ' 5 ) 

Pi 



(PiPb) = iyJ-P b Pi, 
(PaPb) = -vi, 



where we have used the mass-shell condition \pi ± \ = Pi Pi ■ The spinor products fulfill 
the identities (i = Pi,j = pj), 

(ij) = -(ji), 
[ij] = ~ , 



(ij)* = signipiPj) , 

((i + \r\j+)Y = sign(j>° iP °){j + \r\i+), (A.6) 

(ij) [p] = 2 Pi ■ Pj = % , 

(i + \fk\j+) = [ik]{kj), 

(i - = (ik) [k 3 ] . 

B Next-to-leading corrections in the forward-rapidity 
region 

We consider the production of particles p' a , p e , p v in the forward-rapidity region of parton 

Pa, 

yq-Ve^yu^ W, \Pe ± | ^ \Pu ± | ^ \Pq ± | ^ \ P b' ± I , (B.l) 

Momentum conservation ( |A.2|) simply generalizes to, 

Pt - -(Pt+Pt+Pt), (B.2) 
Pb - -Pv ( B -3) 



and accordingly the Mandelstam invariants ( |A.3j ) may be written as 



s 


= 2p a 


■Pb ^ 


(p£ + pt + pt)Pv > 


u 


= 2Pa 


'PV — 


-iPq + Pt + Pt)Pv ' 




= 2p b 


Pk ^ 


-Pk Pv » k = q 7 e 1 p 


4 


= 2p a 


■Pk ^ 


-ipt + Pt + Pt)Pk > 


f 


= 2p b 


Pb' ^ 


-Kl 2 , 



k = q,e,v (B.4) 



30 



to leading accuracy. The spinor products (|A.5|) become 



(PaPb) = -V 7 ! ~ ~yj(j>+ +pf +Pt)Py 

i \ ■ -pi ■ P b '± i \ 

{PaPb> ) = ~l\ — +" Pb' ^ ^ ] APaPb) , 

V Pb' \Pb'J 



/ v . —pi . pi + pi + pi , 

{PaPk} = -i\ —trPk x - -H — + Pk x , k = q,e,v 

V Pk V P k 

(pkPb) = i\j~PbPi - i\jptPv , k = q,e,v 
(pb'Pb) = i^-PbPv - ApvJ > 

/ \ / Pb' / Pk / Pk i 

{PkPb') = Pk ± \ — ~ Pb'\ / Hp - -Pb' W — , k = q,e,is. 
y Pk \ Pb' \ Pb' 
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